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The local persistence probability Pi{t) that a site never becomes active up to time t, and the global 
persistence probability Pg{t) that the deviation of the global density from its mean value p{t) — (p(i)) 
does not change its sign up to time t are studied in a one-dimensional directed percolation process 
by Monte Carlo simulations. At criticality, starting from random initial conditions, both Pi{t) and 
Pg{t) decay algebraically with exponents Oi ft: 9g ^ 1.50(2), which is in contrast to previously 
known cases where 6g < 9i. The exponents are found to be independent of the initial density and 
the microscopic details of the dynamics, suggesting that 8i and 9g are universal exponents. It is 
shown that in the special case of directed-bond percolation. Pi (t) can be related to a certain return 
probability of a directed percolation process with an active source (wet wall). 

PACS numbers: 64.60. Ak, 05.40. +j, 05.70.Ln 



I. INTRODUCTION 

In recent years it has been realized that certain first 
passage quantities in nonequihbrium systems exhibit a 
power law decay with non-trivial exponents |^ ^]. One 
of these quantities is the local persistence probability 
Pi {t) , defined as the probability that a local variable at a 
given space point (normally a spin) has not changed its 
state until time t during a stochastic evolution. In var- 
ious systems it was found that Pi{t) ~ t~^^, where 9i is 
called local persistence exponent. A similar quantity, the 
global persistence probability Pg {t) , defined as the prob- 
ability that the global order parameter does not change 
its sign up to time t is also found to decay as a power 
law with a global persistence exponent 6g. In general the 
exponents 9i and 6g are found to be independent of the 
usual scaling exponents and different from each other. 
The persistence probabilities depend on the history of 
evolution as a whole, and thus it is generally hard to de- 
termine these exponents analytically. Only a few cases 
of exact results are known . 

An important nonequihbrium process which has been 
studied extensively is directed percolation (DP) |11|. A 
large variety of physical systems which undergo a phase 
transition from a fluctuating active phase into an ab- 
sorbing state (i.e. a configuration once reached, the sys- 
tem cannot escape from) belongs to the same universality 
class as DP. Therefore it would be interesting to study 
the persistence probability in DP since this could provide 
some understanding of the non-equilibrium nature of this 
process. 

In this article we present a numerical study of the 
local and the global persistence probabilities in one- 
dimensional (1-d) DP processes. Our results indicate 
that both quantities decay at the critical point as a power 
law with exponents 0; « 0g « 1.5. However, at present it 
is not clear whether these exponents are independent of 
the known scaling exponents of DP. That 9i and 6g are 
the same within numerical errors is surprising since in all 
previously known cases it was found that 9g < 9i M. 



Very little is known about universal properties of per- 
sistence exponents. In a DP process starting from ran- 
dom initial conditions we find the value of 9i to be inde- 
pendent of the initial density, which is in contrast to the 
case of the 1-d Glauber model where 9i depends explic- 
itly on the initial magnetization Q. In addition we find 
these exponents to be the same for different microscopic 
realizations of the DP process, indicating that 9i and 9g 
are indeed universal exponents. 

We also show that the local persistence probability is 
related to certain observables measured in 1-d DP pro- 
cesses with boundaries. Introducing a transfer matrix 
formalism we establish an exact relation between Pi{t) 
and a specific quantity measured in a DP process with 
an absorbing boundary (dry wall). A similar relation is 
found between Pi (t) (which can be seen as a 'first return' 
probability) and the probability R{t) that a DP process 
with a steady active source (wet wall) returns to its ini- 
tial conditions. For a particular realization of DP, called 
bond-directed percolation, this relation can be proven ex- 
actly. Thus the problem of local persistence is related to 
both the dry and the wet wall problems in DP which have 
been discussed in Refs. |l|-|l|. 

The article is organized as follows: In Sect. ||we briefly 
review the DP process. Our numerical results are pre- 
sented in Sect. III. Performing Monte Carlo (MC) sim- 



ulations we measure local and global persistence proba- 
bilities and estimate the corresponding exponents 9i and 
Og. In order to measure 9g more accurately we use a 
block spin method which has been introduced recently 
in Ref. |lj]. In Sect. [V the relation between Pi(t) and 
a specific observable in a DP process with a dry wall is 
proved exactly. The relation between Pi{t) and the re- 
turn probability R{t) in presence of an active source is 
discussed in Sect. 
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II. DIRECTED PERCOLATION 
OVERVIEW 



A BRIEF 



Directed percolation [H] is used as a model for the 
spreading of some, generally non-conserved agent and 
plays a role for certain autocatalytic chemical reactions 
and the spreading of epidemics. In DP models, sites of 
a lattice are either occupied by a particle (active, wet) 
or empty (inactive, dry). In the dynamic processes par- 
ticles can either self-destruct or produce an offspring at 
a neighboring empty site. If the rate for offspring pro- 
duction p is very low, the system will always reach a 
state without particles which is the absorbing state of 
the system. On the other hand, when p exceeds a critical 
value pc another steady state of the system exists on the 
infinite lattice, where the density p{p) of active sites is 
finite. Between the two phases a continuous phase tran- 
sition takes place which is characterized by long range 
correlations. There are various different models for DP, 
e.g. directed site and bond percolation on a lattice [ pl] |, 
cellular automata as the Domany-Kinzel model |16|, and 
the contact process [0 , to name only a few. 

One of the most important properties of directed per- 
colation is its robustness with respect to the microscopic 
dynamics of the system. According to a widely accepted 
conjecture formulated by Janssen and Grassberger [ p^ , 
any transition from a fluctuating active phase into a sin- 
gle, non-fluctuating and non-degenerate absorbing state 
belongs to the DP universality class, provided the dy- 
namical processes are local and characterized by a one- 
component order parameter without special attributes 
like additional symmetries or frozen randomness. The 
DP universality class is characterized by three critical 
exponents, namely the density exponent (i and the scal- 
ing exponents v^^ and v_i. Therefore, if Pi{t) and Pg{t) 
actually decay as a power law, an interesting question 
would be whether 9i and 6g are independent of these 
three exponents and exhibit a similar robustness. 

In the present work we analyze the persistence prob- 
abilities Pi{t) and Pg{t) by Monte Carlo simulations. A 
DP model which is convenient for this purpose is the 
Domany-Kinzel (DK) cellular automaton The DK 
model is defined as follows: a binary variable cr,;(t) — 
0, 1 characterizes the state of site i at discrete time t. 
a — 1 means that the site is active (wet) whereas 
a = means that it is inactive (dry). The automaton 
evolves by a parallel update rule in which the state of 
ai(t + 1) is selected according to transition probabilities 
T{ai{t + l)|(Ti_i(t), (Ti+i(i)). The transition probabilities 
of the DK model are 



r(l|0,0) = 
r(l|0,l) = r(l|l,0)=pi 
r(l|l,l)=P2 
T(0|a-i_i,(Ti+i) = 1 - r(l|(Ti_i,(Ti+i). 



(1) 



10" 



10- 



10 



10"' 



directed bond 

directed site 

Wolfram' s rule 1 i 




10 



10 



10 



10 



t 



FIG. 1. Universal properties of the local persistence ex- 
ponent 9i. The probability Pi{t) that a site never becomes 
active up to time t is measured in a Domany Kinzel model 
with 1024 sites averaged over 10^ independent runs. Starting 
with the initial density p(0) = 0.8 various critical points are 
examined. The results indicate a power law decay which is 
consistent with 9i — 1.49(2) in all cases. The bold line marks 
slope -3/2. 



where a continuous phase transition takes place. It is 
assumed that for p2 < 1 the critical behavior along the 
whole line Pi,c(P2) is the same as that of DP. In partic- 
ular, the phase transition line includes three particular 
points, namely | |T9|] : 



1. directed bond percolation: pi = 0.644701(1), 

P2 = 2pi 



Thus the model is controlled by two parameters pi 
and p2. For fixed P2 < 1 there is a critical value pi,c 



2. directed site percolation: pi = P2 = 0.705485(5) 

3. Wolfram's rule 18: pi = 0.8092(1), p2 = 

As a necessary condition, a universal critical exponent 
should not depend on the choice of the transition point 
used in a simulation. 

Another special case is compact directed percolation 
where pi — 1/2 and p2 = 1- Here the dynamics of the 
model belong to a different universality class, namely 
that of annihilating random walks or, equivalently, the 
Glauber Ising model. In fact, it turns out that in the 
case of compact DP the results of Rcfs. arc recov- 

ered. 



III. MONTE CARLO SIMULATIONS 

A. Local persistence probability 

The local persistence probability Pi (t) is defined as the 
probability that in a DP process starting from random 
initial conditions a given inactive site does not become 
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FIG. 2. Universal properties of tlie local persistence expo- 
nent 9i. Piit) is measured for directed bond percolation with 
various initial densities. For p(0) > 0.4 an algebraic decay 
with 9i — 1.50(2) is observed. In the case of very low initial 
density there seems to be a longer transient until the final 
power law decay is reached. The bold line indicates the slope 
-3/2. 

active up to time t. This quantity is nontrivial because 
it can be seen as an infinite-point correlation function in 
the direction of time. 

Notice that in contrast to Ising systems we define Pi (t) 
as the probability for a site not to become active rather 
than not to flip. This is because there is no symmetry 
between active and inactive sites in DP: Since active sites 
can spontaneously turn into inactive sites, the probability 
for a site to remain active up to time t decays exponen- 
tially. On the other hand P; (t) , defined as the probability 
for a site to remain inactive, is expected to decay more 
slowly since inactive sites can only become active in pres- 
ence of an active neighboring site. 

The qualitative behavior of Pi (t) depends on the per- 
colation probability. In the inactive phase the density 
of active sites decays exponentially fast until the system 
enters the absorbing state. Therefore a finite fraction 
of sites remain inactive so that Pi{t) saturates at some 
constant value. In the active phase an infinite percolat- 
ing cluster emerges. Since the size of inactive islands in 
this cluster is finite, Pi{t) will decay exponentially in that 
case. At the critical point, however, there is no charac- 
teristic length scale, wherefore we expect Pi{t) to decay 
as a power law. 

We performed Monte Carlo simulations for various 
transition points and different initial densities. In all 
cases a power law decay is observed. The local persis- 
tence exponent 9i « 1.50(3) is found to be independent 
of the choice of the transition point as well as the initial 
density (see Figs. |l|and||). Similar results (not reported 
here) are obtained in a contact process [E7| which is a 
model for DP with random sequential dynamics. There- 
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FIG. 3. Finite-size scaling function for the local persis- 
tence in the absorbing phase. The scaled persistence prob- 
ability Pi{t)t''' is measured at criticality for system sizes 
L = 8, 16, 32, . . . , 1024 and agains scaled time t/L^ . The best 
data collapse is obtained for 9i — 1.50(1). 

fore our findings suggest that 9i is indeed a universal 
exponent in DP. 

In addition we examined the scaling properties of the 
local persistence probability. According to the usual scal- 
ing theory of DP [|ll| , we expect Pi (t) to scale like 

Piit,L,e)^t-'' fie''H,L-H), (2) 

where e — \p~Pc\ measures the distance from the critical 
point, z = v\\lv]_ is the dynamical exponent, and / is a 
universal scaling function with the asymptotic behavior 

/(0,0) = const 

f{x,0) - x'^' for x-^oo (3) 
f{0,y) ^ y^' for y ^ oo. 

We verified the finite-size scaling of Eq. (p|) at critical- 
ity. In Fig. ^ the scaling function /(O, L'^t) is shown for 
various system sizes. The best data collapse is obtained 
for 9 = 1.50(1) which is our the most precise estimate 
for the local persistence exponent. Attempts to relate 9i 
to the DP exponents /3, i'±, have failed wherefore we 
believe that 9i is an independent critical exponent. 

B. Global persistence probability 

The global persistence probabihty Pg{t) is usually de- 
fined as the probability that the global order parameter 
(analogous to the total magnetization in Ising Model) 
does not change its sign up to time t. It has been studied 
in various models as, for example, in the Ising model 0| 
and models in the parity conserving class 0]. In all 
cases it was found that Pg{t) decays algebraically with 
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an exponent 9g independent of 9i and the other scaling 
exponents of the system. 

In DP, however, the global order parameter - the den- 
sity of active sites p{t) ~ is a non-negative quantity and 
therefore the above definition is not applicable. Instead 
we may consider the probability that the deviation of the 
order parameter from its mean value /^p{t) = p{t) — {p{t)) 
does not change its sign up to time t. As in the case of lo- 
cal persistence it turns out that this probability depends 
on the sign of Ap(t) since there is no symmetry between 
active and inactive sites. In fact, the probability for the 
deviation to remain positive up to time t decays expo- 
nentially fast whereas the corresponding probability for 
negative deviation decays more slowly. We will there- 
fore define Pg{t) as the probability that ^p{t) < up to 
time t. 

Usually it is very difficult to measure 9g by Monte 
Carlo simulations, the reason being that each run pro- 
duces only one data point (the first passage time) whereas 
the measurement of 9i produces a number of data points 
of the order of the system size L. Hence a large number 
of runs is required to obtain good statistics. However, 
it has been shown recently that this problem can be 
circumvented by introducing spin-block persistence prob- 
abilities Pm{t) which are defined as follows: Consider 
blocks of m sites and define a block density pm (t) as the 
average fraction of active sites in each block. Pm{t) is 
defined as the probability that Apm{t) = pm{t) — {p{t)) 
remains negative up to time t. Obviously the spin-block 
persistence probability Pm{t) connects the special cases 
of local and global persistence, namely 

P,{t)^Pl{t), P^{t)=Pg{t). (4) 

It was observed in Ref. that Pm{t) in a Glauber 
model first decays as Pg{t) ~ t^^^ and then crosses over 
to the power law decay Pm{t) ~ Pi{i) ~ t~^^ , where the 
crossover time grows with the box size. Since the number 
of such boxes is of the order L/m, this method allows to 
measure 9g much more accurately. 

We have measured the Pm{t) for block sizes m — 
1,2,4,8,16. As can be seen in Fig. ^, P„i{t) decays as 
t~^"^ where the exponents 9m — 1.50(4) seem to be in- 
dependent of m, suggesting that the global and the lo- 
cal persistence exponents are identical. This is surpris- 
ing since in all previously known cases 9g was found to 
be smaller than 9i. For example, in the Id Ising model 
0,, = 1/4 and ^,=3/8 

In order to support this observation, we also measured 
the global persistence probability Pg (t) over two decades 
in time (see bold line in Fig. ^). Although this measure- 
ment is not very accurate, there is no indication for 9g 
to be smaller than 9i which supports the previous result. 
At present we have no explanation why the exponents 9i 
and 9g appear to be equal in DP. 
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FIG. 4. Global persistence probability and Block spin 
method: Decay of the persistence probability Pm (t) for block 
sizes m = 1,2,4,8,16. From the slopes we obtain the esti- 
mates Oi = 1.52(2), 02 = 1.48(2), 94 = 1.49(2), Os = 1.52(2), 
9ie = 1.49(2). The bold line represents a direct measurement 
of Pg{t) (shifted vertically by a factor of 10). The dashed line 
indicates slope -3/2. 

C. Power law versus stretched exponential 

In certain reaction diffusion models it was shown that 
the Pi{t) decays as a stretched exponential function 
rather than a power law Q . In numerical simulations it is 
sometimes difficult to distinguish between stretched ex- 
ponential and power law decay. In order to verify whether 
Pi (t) truly decays as a power law one can use a heuristic 
argument as indirect test |^ : For a site which was never 
active up to time t to become active in the next time step 
t + 1, it is necessary that a neighboring site is active at 
time t. Hence Pi{t) changes according to 

lp,{t) = -Pi{t)psit) (5) 

where Ps{t) is the probability for finding an active site 
near a site which was never active until time t. By in- 
tegrating Eq. (^ one can easily see that Pi{t) decays as 
a power law only if Ps{t) t^^. On the other hand, if 
Ps{t) decayed as t~" with a 1, Pi{t) would decay as a 
stretched exponential. Thus by measuring the exponent 
a in a Monte Carlo simulation we can verify the observed 
power law decay for Pi{t). Our results (see Fig. ||) are 
consistent with Ps{t) ~ t^^, supporting that Pi(t) actu- 
ally decays algebraically. 



IV. RELATION TO DIRECTED PERCOLATION 
WITH AN ABSORBING BOUNDARY 

In this section we prove that the local persistence prob- 
ability Pi{t) is exactly equal to the expectation value of 
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FIG. 5. Verification of the power law decay of Pi{t). The 
probability Ps{t) for finding an active site near a site which 
was never active until time t is plotted as a function of time. 
It is observed that ps(t) ~ t"" with a = 1.01(2), indicating 
that Pi{t) decays indeed algebraically (see text). 



a specific observable in a DP process with an absorbing 
boundary. A boundary is called absorbing if all bonds 
across the boundary are cut. In a 1-d system such an 
absorbing boundary can be introduced by forcing a par- 
ticular site to be inactive during the whole time evolu- 
tion, i.e. an absorbing boundary can be understood as a 
dry wall. The effect of an absorbing boundary in a 1-d 
DP process has been studied recently in Refs. [p"2|-p^, 
and it is therefore interesting to investigate the relation 
between the two problems. 

Let us consider a DK model with L sites and periodic 
boundary conditions. Denoting by \a) = {cti, (T2, . . . , (Tl} 
basis vectors in configuration space, the transfer matrix 
T of the DK model is defined by 

L 

(a'(t + l)|T|a(i)) =[]r((7,(t + l)|a,_i(i),a,+i(i)) . 



(6) 



Furthermore let us define vectors for the absorbing con- 
figuration |0), the state \a) where all sites are active, and 
the sum over all configurations |1): 



|0) = {0,0, 
|a) = {l,l, 

|1) = 



(Ti,0-2----,(7l 



0,0} 
1,1} 



(7) 



Using this notation the transfer matrix obeys 
T|0) = |0), (1|T=(1|, 



(8) 



i.e. the absorbing state is a ground state of the DK model 
and the transfer matrix conserves probability. Let us now 



select an arbitrary site j and define local operators So, 
Po, and W by 



(a' I Sola) 
(^'|Pok) 



L 

i=l 



(9) 



(a'|W|a)=T(0|a,_i,a,+i)n(5,-,., . 

i=l 

So turns site j into the inactive state, Po projects onto 
states where site j is inactive, and W is a diagonal weight 
operator to be explained below. Notice that So conserves 
probability whereas Po does not. Using these notations 
we can now express the local persistence probability Pi (t) 
as 



PKt) = (l|(PoT)*N), 



(10) 



Here \i) denotes the initial probability distribution where 
the average over different initial configurations is taken 
into account. For example, the average over random ini- 
tial conditions with particle density p corresponds to the 
initial distribution 

L 

(71,(72, ...,CTi i—l 

In Eq. (1^ the projector Po removes all space-time his- 
tories which would turn site j into the active state giving 
just the local persistence probability Pi{t). Notice that 
this expression is already properly normalized. 

Now consider a DP process with an absorbing bound- 
ary inserted at site j. This process is described by a 
different transfer matrix T where site j is forced to be 
inactive: 



SoT. 



(12) 



One can easily show that T is related to the full transfer 
matrix T byPoT = TW which implies that 



P,{t) ^ {1\{'TWY\^) . 



(13) 



Therefore the local persistence probability Pi{t) is ex- 
actly equal to the expectation value of the diagonal op- 
erator W measured before each update in a DP process 
with an absorbing boundary. Eq. (n3[) can be written as 



Pi{t)^{l[T{0\a,.,{t'),a,+^{t')) 



(14) 



where () denotes the average over many independent re- 
alizations of a DP process with an absorbing boundary 
at site j combined with an independent average over ini- 
tial configurations according to the probability distribu- 
tion Note that T{0\(jj^i(t'),aj+i(t')) refers to the 
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transition probability of the DK model without dry 
wall at site j. 

Because of the dry wall at site j, aj^i{t) and (Tj+i(t) 
are uncorrelated in an infinite system. Assuming that 
events where aj-i{t) = (Jj+iit) = 1 are rare we may 



line A 



therefore approximate Eq. (|14[ ) by 



exp 



t-i 

t'=0 



(15) 



where 7 = ~21og(l —pi). Thus Pi{t) is related to the 
integrated activity next to the boundary. As mentioned 
before, the surface activity of a 1-d DP process with an 
absorbing boundary has been carefully analyzed by series 
expansions and MC simulations in Refs. It was 

observed that the activity next to the boundary decays 
according to a power law 



(16) 



Furthermore, it has been conjectured that /?' = — 1 
which means that 
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FIG. 6. Mapping of the local persistence problem onto 
a directed bond percolation process with a steady source in 
the center. The figure shows a particular realization of open 
bonds (solid lines) on a tilted square lattice (dotted lines) in 
directed bond percolation. Active sites are marked by bold 
dots. Initially all sites are in the active state (line A). Activity 
then percolates through the system without touching site j 
(line B) . In reverse time direction this means that there is no 
open path backwards from line B to line A. 



(17) 



\t'=0 



If the average () in Eq. (15) commuted with the expo- 
nential function, Eq. ( [l7| ) would imply that Pi{t) decays 
asymptotically as a stretched exponential 



Pi{t) ^eM-lt'^"^^) 



(18) 



rather than a power law. However, our numerical results 
strongly sugge st tha t Pi (t) does indeed decay as a power 
law (see Sect. IIIC ). In fact, a stretched exponential of 
the form ( |l8| ) is in obvious contradiction with the simu- 
lation data. Thus the average operation () certainly not 
commutes with the exponential function. Nevertheless 
we are left with a puzzle: For both Eqs. (|l^) and ( |l7| ) 
to decay algebraically, a delicate mechanism in the ex- 
ponential function is needed, i.e. higher cumulants of 
the integrated surface activity have to match in a very 
specific way. 



V. RELATION TO A DIRECTED PERCOLATION 
PROCESS WITH AN ACTIVE SOURCE 

We will now show that the local persistence probability 
Pi{t) can be related to a return probability in a Id DP 
process with a pointlike active source. More precisely 
we will show that Pi{t) is equal to the probability for a 
DP process with an active source to return to a state 
where all sites except for the source are inactive. We will 
prove this relation exactly in the case of directed bond 
percolation. 

The mapping can qualitatively be understood as fol- 
lows: Consider a particular realization of a DP process 



(a single MC run) starting from initial conditions where 
all sites except for site j are active. Let us assume that in 
this particular realization site j has never become active 
up to time ti. Such a realization is shown schematically 
in Fig. ^ Obviously there are no open paths from the 
horizontal line t — (line A in Fig. ^) to the vertical line 
at site j given by < ? < ti (line B). Conversely, there 
is no open path backwards in time from line B to line A. 
In the special case of directed bond percolation we may 
now consider a DP process in reverse time direction using 
the same realization of open and closed bonds. Further- 
more, let us assume that we force site j to be active along 
line B, i.e. we impose an active source at this location. 
If there is no open path from B to A, activity will not 
percolate from line B to line A in the 'time-reversed' pro- 
cess. In other words, the reverse process returns to its 
initial condition where all sites except for the source are 
inactive. 

We now prove this mapping for the special case of di- 
rected bond percolation. More precisely we show that 
in this case the probability R{t) of a DP process with 
an active source to return to its initial condition where 
all sites except for site j are inactive is exactly equal to 
the persistence probability Pi{t). Using the notations of 
Sect. let us introduce two further operators Si and U: 



L 

(a'lUlfi) =[](!- -^.^i^..,!). 



(19) 



i=l 



The operator Si turns site j into the active state whereas 
the operator U is a symmetric transformation matrix. 
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One can easily verify that for p2 = Pi(2 — pi), i.e., for 
directed bond percolation, the following relation holds: 



UT = U . 

Furthermore we have 



UPo = S(U, 



USo 



sfu. 



(20) 



(21) 



By commuting the matrix U to the right and transposing 
the resulting expression, we can now rewrite the local 



persistence probability Pi{t) in Eq. (13) as 



m) = {l\{PoTy\^) 

= (0|U(PoT)*So|a) 
= (0|(SfT^)*S;^U|a) 
= (a|USo(TSi)*|0) 
= (0|So(TSi)*|0)=i?(i), 



(22) 



where we assumed the initial condition \i) ~ So|a) 
for the persistence measurement in which all sites ex- 
cept for site j are active. We also used the relations 
(1| = (0|U and (a|U = (0|. The resulting expression 
(0|So(TSi)*|0) = R{t) is precisely the return probability 
of a DP process to its initial condition with an active 
source at site j, which completes the proof. Notice that 
similar arguments were used in Ref. in order to derive 
the persistence exponent in the Glauber model. 

Although this proof holds only for the case of directed 
bond percolation, the relation seems to be more general. 
In fact, we verified numerically that Pi{t) and R{t) ex- 
hibit the same power law behavior for various transition 
points in the DK model. 



VI. CONCLUSIONS 

In the present work the problem of local and global 
persistence in directed percolation has been studied nu- 
merically. The results suggest that both the local and 
global persistence probabilities exhibit an algebraic de- 
cay in time at the critical point. The corresponding crit- 
ical exponents 9i and 6g are universal and independent of 
the initial density. The universality of these exponents is 
stronger than in the Ising model where an explicit depen- 
dence on the initial density is found for 9i. The reason 
might be an exponentially fast decaying memory of ini- 
tial conditions in DP since active sites can spontaneously 
become inactive. Surprisingly, both 9i and 9g seem to 
be equal within numerical errors which contrasts to pre- 
viously known cases where it was found that 9g < 9i. 
We carefully analyzed our data in order to rule out any 
stretched exponential decay of the persistence probabil- 
ity. 

In order to consider the problem in a more general 
context, we related the local persistence probability to 
certain observables in a DP process with a dry and a 



wet wall, respectively. Introducing a transfer matrix for- 
malism these relations were proven exactly for particular 
realizations of DP. 

Various questions remain open. First of all, we do not 
know whether 9i and 9g are in fact independent from 
the usual critical exponents in DP. The numerical value 
9i « 1.50(1) suggests that the exact value could be 3/2, a 
possibility which cannot be ruled out as an integer expo- 
nent was also observed in a DP process with an absorb- 
ing boundary |]l^ , |l3| . Furthermore, it is not clear why 
9i and 9g appear to be identical in DP. Finally, it would 
be interesting to investigate the same problem in higher 
dimensions. 
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